Charge separation in Reheating after Cosmological Inflation 
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New aspects of parametrically resonant heating of a relativistic scalar 0(2)-symmetric self-interacting field 
are presented. This process is a candidate for reheating at the end of the early-universe epoch of inflation. 
Although a model with a fully symmetric ground state is used, transient, metastable spontaneous symmetry 
breaking can be observed. This manifests itself in the form of persistent regimes of opposite and, inside these, 
uniform charge overdensities separated by thin lines and walls similar to topological defects, in two and three 
spatial dimensions, respectively. The configuration is found to correspond to a non-thermal fixed point of the 
underlying dynamic equations for correlation functions which prevents thermalisation over an extended period 
of time. Our results establish a link between wave-turbulent phenomena and the appearance of quasi-topological 
defects in inflaton dynamics. 
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In cosmological models of the universe, reheating describes 
the epoch starting at the end of inflation Q. During this epoch 
the potential energy of the inflaton field is redistributed into a 
homogeneous and isotropic hot plasma of particle excitations. 
These become a substantial part of the further expanding uni- 
verse. Simple models describing reheating after inflation in- 
voke self-interacting scalar fields. One of the popular sce- 
narios involves the parametrically resonant amplification of 
quantum fluctuations of the macroscopically oscillating infla- 
ton field. The amplified modes represent the emerging matter 
content of the universe Various theoretical approaches 

have been proposed to model reheating. As both, the infla- 
ton and the amplified modes are strongly populated, classi- 
cal field simulations can be applied to describe their evolution 
f4-6|. Alternatively, Kadanoflf-Baym dynamic equations, as 
derived from 2PI eflTective actions in nonperturbative approx- 
imation, can describe the resonant excitation and the ensuing 
thermalisation |7 1. The exponentially fast excitation process is 
followed by a slower equilibration, possibly with transient tur- 
bulent behaviour transporting the energy deposited in the low- 
momentum modes of the system to higher momenta. Classi- 
cal field simulations and scaling solutions of kinetic equations 
were used to analyse possible turbulent evolution during re- 
heating and thermalisation . 

Recent analytical predictions of new nonthermal fixed 
points of the field dynamics |10| provided an extension of 
the above turbulence scenarios into the regime beyond per- 
turbative quantum-Boltzmann approximations (HI [121. Be- 
ing confirmed by numerical simulations (TOl [131 these phe- 
nomena have triggered work in diflTerent areas (121 [HI [TSl , 
in particular as they have the potential to strongly alter relax- 
ation time scales. Here we demonstrate that these nonthermal 
fixed points are characterised by metastable configurations of 
the underlying inflaton field with a soliton-like charge distri- 
bution. We study parametric resonance in a globally 0(2) or, 
equivalently, U(l) symmetric relativistic scalar field theory. 
Shortly after the resonant excitations have set in we find spa- 
tial separation of charges. Both, charge and anticharge over- 



densities become uniformly distributed within slowly vary- 
ing regions which are separated by sharp boundary walls of 
grossly invariant thickness. These walls have a character sim- 
ilar to topological defects and appear for generic initial condi- 
tions. Metastable, non-topological defects have been studied, 
mostly in the context of non-linear single- scalar field theo- 
ries, with or without a broken- symmetry ground state, where 
they are commonly termed oscillons |[T6l . In the present paper 
we connect quasi-stationary soliton-like charge separation in 
a two-component, 0(2)-symmetric model to non-thermal sta- 
tionary scaling solutions of quantum field dynamic equations 
which are known to exist for a far more general class of mod- 
els. In this way a link is established between wave turbulence 
phenomena as discussed, e.g., in ll8UT4l and long-lived quasi- 
topological structures in the inflaton field. The very general 
concept of critical phenomena far from thermal equilibrium 
can provide a new way to classify nontrivial metastable field 
configurations. 

The action of the model considered here is given by 

5= J Ajf|if(5,^)2- (5,^)2 -mVl-^(^¥} (1) 

where (f^ = (fia(x)(fa(x), {di(pf = di(fa(x)di(fa(x), sums over 
/ = 1, . . . , J and a = 1, . . . ,N, N = 2, implied. To describe the 
dynamics of the scalar field (p in the rapidly expanding uni- 
verse one works in conformal coordinates in which the metric 
is defined through ds^ = a(t)^(dt^ - dx^), with a scale param- 
eter a(t) depending on the conformal time t. We choose the 
bare mass scale to vanish, = 0, such that the equation of 
motion for the reseated field (p' = a(f reads (in d = 3): 



During the parametric reheating epoch the universe is close to 
being radiation dominated |4|, such that the scale dependent 
term can be neglected, a = 0. Hence, the equation of motion 
for the reseated field O = (f' /(f^, with (f^ = (f\t = 0) the initial 
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Figure 1 : (Color online) Occupation-number spectrum ^ of the system in J = 2 (left panel) and d = 3 (right) spatial dimensions at different 
times t (in lattice units). Averages were taken over 69 (84) runs on a 960^ (256^) grid. The double-log scale exhibits the bimodal power-law 
n(k, t) ~ k~'^. For d -IwQ find /c ^ J + 1.5 in the infrared, and thermal scaling, k ^ 1, for large k. For J = 3 we also observe, at late times and 
intermediate momenta, the UV weak wave turbulence exponent k ^ 1.5. Note that the physical particle number spectrum is 6Nn(k, t)/A. See 
the main text for further details. 



value of the inflaton field, can be approximated as 

[^^-A + O^jo, = (3) 

where the coupling has been absorbed by the rescaling 
yJXlWip'^x X and ^|X|6N^p'^t t. Due to the global 0(2) 
symmetry of the action the charge-current density 

Ux) = Oi(x)^^02(x) - ^2(x)d^^i(x). (4) 

is conserved, = 0. 

To induce parametrically resonant reheating, we start our 
simulation from a configuration where only the zero mode of 
the inflaton field is populated. Fluctuating non-zero momen- 
tum modes act as seeds for the ensuing instabilities. Choos- 
ing = and A > 0, which corresponds to an equilib- 
rium configuration in the symmetric phase, subsequent oscil- 
lation of the inflaton field induces parametrically resonant ex- 
ponential growth of certain modes. Scattering between these 
modes causes the entire spectrum to fill up. Our simula- 
tions of Eq. ([3| were performed on a cubic space-time lattice 
of gridpoints using leap-frog discretisation and periodical 
boundary conditions in d = 2 and 3 spatial dimensions. The 
time evolution of the system after the instabilities have set in 
was found to be best described by a turbulent stage followed 
by thermalisation on the largest time scales. The turbulent 
phases feature universal scaling of the (solid-)angle-averaged 
momentum-space occupation numbers 

n(kj)=^ / ^^"'^kV<l^.^(k,0P><|O(k,0P). (5) 

From standard analytical studies based on quantum Boltz- 
mann equations ifTTl . one expects to find, during these phases, 
the above radial occupation number to scale as n(k, t) ~ k~'^ 
where the specific power-law exponent k reflects weak wave 
turbulence. The distribution is approximately time indepen- 
dent, n{k, t) = n(k), for k within the so-called inertial range. 



but in contrast to thermal equilibrium, energy is flowing 
through this range, between sources and sinks, establishing 
a non-equilibrium stationary state. Diflferent exponents k are 
possible, depending on whether the radial flow of energy is in- 
dependent of k, P(k) = P, or the radial flow of (quasi- )particle 
number, Q(k) = Q L17.1 . In summary, one finds the possible 
exponents ifTUlfHlfTTl 
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Kp = d-\- -(z-2), KQ = d-\-z-3. (6) 

where z is the homogeneity index of the dispersion co(sk) = 
s^aj(k) of mode k. Note that in thermal equilibrium one ex- 
pects Kth = z, corresponding to the Rayleigh- Jeans law. Given 
positive K, occupation numbers grow large in the infrared (IR) 
and the Boltzmann perturbative approximation breaks down. 
In |[T0l[T2l, an expansion of the two-particle irreducible (2PI) 
eflfective action to next-to-leading order in the inverse 1/N 
of the number of field components was used to extend the 
wave-turbulence analysis into the IR regime. In this scheme, 
bubble-chain diagrams contributing to the self-energy in the 
dynamic Dyson equation are resummed, resulting in an eflfec- 
tively renormalised many-body coupling. With this, scaling 
exponents can also be derived in the IR regime: 

4^ = J + 2z, Kf = d^z (7) 

As was first brought forward in 1 10] for the case of an 0(4) 
model, numerical simulations of parametric resonance con- 
firm both, the IR and the UV scaling exponents. Correspond- 
ing results for the model studied here are shown in Fig. [T] for 
d = 2 and d = 3 at different instances of time t. 

As was shown in 1 17], for our model ([T]), in wave turbu- 
lence, either energy flows with rate P towards larger k or 
quasiparticles flow with rate Q towards lower k. In our simu- 
lations, the resonantly induced momentum distribution is ob- 
served to act, during the further evolution, as a source at inter- 
mediate momenta, k ^ 0.1.. 0.3, see Fig.[T] Our model implies 
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Figure 2: Modulus |O(x,0l = [^i(x,0 + O^(x,0]^^^ at lattice time 
? = 9 X 10^. Shown is its distribution over the two-dimensional 960^ 
lattice used in our simulations. "Defect" lines with substantially re- 
duced p = jo (cf. (|4|), contrast with a fairly constant bulk pbuik - 0.7. 
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Figure 3: Charge density p(x,t) = 7o(x, (cf. (|4])), for the same 
configuration as shown in Fig. [2] Note the evenly distributed charge 
overdensities with opposite sign, separated by lines which include 
the "defect" lines shown in Fig.|2] 



z = 1 such that we expect from (|7]), in the infrared, the expo- 
nent a:^ = J + 1, assuming a constant quasiparticle number 
flow towards the IR. We note that in contrast to cases with 
smaller l/N discussed in iW.TTl we find, in the IR, k ^ 3.5 
(d = 2) and 4.5 (d = 3). In the UV, assuming a constant flow 
P from intermediate to larger momenta one expects, from ([6]), 



/Cp^ = d -3/2. Numerically we find that in the UV the occu- 



pation number distribution approaches, at large times, thermal 
scaling Kth = 1 . For J = 3 we also observe, at late times and 
intermediate momenta, the expected UV weak wave turbu- 
lence exponent k ^ 1.5. 

Having recovered scaling behaviour we present, in the fol- 
lowing, our central result: Looking at the real-space struc- 
ture of the emerging critical configuration we find patterns 
similar to topological defects giving rise to quasi- stationary 
charge separation. We find a correspondence between the ap- 
pearance of the strong scaling exponent and of the defect- 
separated charge patterns. In Fig. [2] we depict, for J = 2, a 
typical real-space configuration in the turbulent stage, plot- 
ting the modulus of the 0(2) scalar field, |0(x, 01 = [<I>i (x, t) + 
(^l(x,t)]^^^. Localised regions appear, specifically "defect" 
lines where the absolute value of the field is much smaller 
than its average. Fig. |3] shows the corresponding charge den- 
sity p(x, t) = jo(x, t), Eq. (|4]). Clearly, both uniform charge 
and anti-charge overdensities appear within distinctly sepa- 
rated regions, showing only small fluctuations as compared to 
their bulk values. This separation of charges is confirmed by 
the histogram of local charge densities shown in Fig. [4] We 
found similar histograms for the case of the 0(N) model ([T]) 
with N = 3,4 which we will study in more elsewhere. 

Moreover, we see that the "defect" lines in Fig. |2] are ly- 
ing along the boundary limits between the regions of opposite 
charge. Note that only part of the boundaries is clearly visible 
as "defect" lines. Following the time evolution of the defect 
lines © I we find that they appear and disappear with a pe- 
riod of r/2 ^ 6 determined by the eflTective mass meff = 2n/T 
of the mean field theory, at every point along the boundaries 
between the oppositely charged domains. The explanation for 
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Figure 4: Histogram of the charge density distribution of the system 
in J = 2, at points of time for which Fig. [T] shows the corresponding 
spectra. We find that a clear charge separation goes together with the 
strong IR scaling with /c = 3.5. 



the appearance of these lines is as follows: The charge den- 
sity p can be interpreted as field angular momentum in the 
N = 2 dimensional field space. Hence, J-dimensional con- 
figuration space can be divided into regions according to the 
sign of p. On the lines separating these regions the local field 
variable has neighbours of opposite charge which therefore 
possess opposite "circular polarisation" in field space. Thus, 
at a particular point on the boundary the rotations add to a 
linear oscillation of the field through zero, giving rise to spa- 
tially and temporally oscillating defect lines. While such an 
evolution could similarly occur in a free theory (A = 0) in- 
teractions are responsible for rendering the transition between 
the opposite charges sharp as can be seen in Fig. [2] 

Following the evolution of the system starting from the ho- 
mogeneous initial configuration we find the defect lines to 
appear once the instability has saturated. The lines can re- 
connect with each other while they propagate slowly through 
configuration space, such that the number of spatially con- 
nected charged regions changes in time. On the average, the 
total length of the defect lines as well as the number of con- 
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Figure 5: Charge density p(x, = 7o(x, (cf. (|4])) for the same run 
as shown in Fig. [5] at the times t - 300 (upper left), ^ = 10^ (upper 
right), t = 10^ (lower left), ^ = 6 X 10"^ (lower right). The scale of the 
largest structures is found not to depend on the total lattice size. 



nected regions is found to decrease with time, see Fig.|5] Dur- 
ing the late stage of the evolution the separation of opposite 
charges gradually disappears while the system finally ther- 
malises. This can also be seen in the histogram in Fig. |4] 
where the central well gradually disappears while the spec- 
trum n{k, t) approaches thermal equilibrium, see Fig. [T] We 
note that the sharp step between the opposite charges remains 
(Fig. [Slbottom right panel). This, as well as the histogram 
in Fig7|4] suggest that topological defects occur according to 
the homotopy group of the charge density space C instead of 
the order parameter space. Although C, for the 0(2) model, 
equals the real numbers, our results show, that the system re- 
mains, for a long time, in a subspace which is topologically 
equivalent to 5*0 = {-1, 1} ^ C. Hence, one expects defect 
lines (walls) to separate regions corresponding to the two un- 
connected parts of C. 

For J = 3, starting from an equivalent initial configuration, 
we find bubbles of opposite charge separated by thin walls 
(Fig.|6]) [18] as well as charge density histograms similar to the 
case d = 2. These bubbles are seen to appear in coincidence 
with the strong IR scaling in the spectra (Fig.[T] right panel). 

We emphasise that the structures found are distinctly dif- 
ferent from the known topological defects appearing in 0(2) 
theories such as vortices ind = 2 and vortex strings ind = 3 
which one obtains in the case of equilibrium symmetry break- 
ing, i.e., < lO. The bubble walls reported here form for 
vanishing or positive mass squared. 

Considering other potential cosmological consequences of 
our findings, notice that, in a diff'erent context, the U(l)- 
symmetric scalar field theory can also be understood as an in- 
gredient of Affleck-Dine baryogenesis 1 19 |, where the scalar 
field carries a baryon charge and is thought to be a supersym- 
metric partner of standard-model fields. We emphasize that, 
in contrast to such models, there is no symmetry breaking im- 
plemented in our model from the outset. The breaking rather 
occurs as a dynamical, transient but quasistationary eff'ect. 

To summarize, we have studied parametric reheating in an 




Figure 6: Bubbles of approximately constant and opposite charge 
separated by thin walls on a 64^ lattice at ? = 1500. Plotted are 
surface lines which are fitted to lattice points having |Op < 0.05. 



0(2) symmetric scalar field theory in two and three spatial di- 
mensions. We have found structures similar to (quasi-) topo- 
logical defects to emerge on intermediate time scales, break- 
ing charge symmetry locally while the classical action im- 
plies an equilibrium configuration without symmetry break- 
ing. These structures were found to appear in direct corre- 
spondence with strong wave-turbulent scaling behaviour of 
occupation numbers in the IR regime of low momenta. This 
establishes a strong link between wave-turbulence and quasi- 
stationary defect-like pattern formation of inflaton dynamics 
in scenarios of reheating. It is an interesting question whether 
similar structures exist in more complicated theories, such as 
gauge theories. 
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